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ABSTRACT 

Continuous  radio-wave  telemetry  is  required  during  planned  tests  of  direeted-energy 
weapons  systems  in  order  to  characterize  in  situ  the  effects  of  laser  irradiation  on  different 
target  materials.  Unfortunately,  the  incident  radiation  ean  cause  disruption  of  the  radio  signal 
during  the  direeted-energy  testing.  Several  phenomena  associated  with  direeted-energy 
impact  can  lead  to  communication  path  losses,  such  as  ablation,  charged  particle  emission, 
charring,  and  chemical  changes  in  the  target  materials.  Direeted-energy  impact  on  the  target 
material  leads  to  target  heating  and  consequent  ablation.  A  numerical  model  has  been 
developed  to  describe  the  laser  induced  ablation  of  metal  surfaces.  The  model  describes  the 
absorption  of  the  laser  energy  by  the  metal  and  the  resulting  temperature  rise  in  the  surfaee. 
This  temperature  rise  then  induees  ablation  of  the  target  material.  Results  for  an  aluminum 
target  irradiated  with  a  KrF  laser  were  obtained.  Temperature  profiles  in  the  target  material 
and  surfaee  temperature  ehanges  are  presented  along  with  the  ablation  rate  as  a  function  of 
time  as  the  aluminum  target  is  irradiated.  This  report  presents  results  for  eases  when  laser 
energy  absorption  by  the  plasma  plume  created  above  the  surfaee  is  not  significant. 


INTRODUCTION 

Direeted-energy  impact  on  the  target  material  leads  to  target  heating  and  consequent  ablation. 
The  ablated  material  continues  to  expand,  compresses  the  air  and  forms  a  shock  wave  in  the 
air1.  Within  the  plasma  and  near  the  shock  front,  the  combination  of  high  air  temperatures 
and  the  ablated  material  causes  complex  chemical  reactions  leading  to  ionization, 
recombination,  and  excitation.  In  addition,  chemical  reactions  between  the  ablated  material 
and  air  molecules  are  possible  at  these  elevated  air  temperatures. 

There  have  been  a  number  of  numerical  models  proposed  to  analyze  the  laser-solid 
interaction  and  ablation  process111"14'".  We  developed  a  model  based  on  a  kinetic 
description  of  the  Knudsen  layer  and  a  hydrodynamic  description  of  the  collision-dominated 
plasma  region.  Preliminary  analysis  of  the  ablation  rate  of  various  targets  subject  to  direeted- 
energy  impact  was  performed.  The  ablation  rate  depends  on  the  surfaee  temperature  as  well 
as  plasma  density  in  the  target  vicinity^.  In  this  present  paper,  we  report  on  an  initial  attempt 


to  apply  this  ablation  model  to  a  general  ease  of  an  aluminum  target  being  irradiated  by  a 
KrF  laser  in  vacuum  conditions.  The  model  prediets  the  ablation  rate  and  the  depth  of  the 
resulting  crater,  as  well  as  the  temperature  profile  in  the  target  material. 


MODEL 

A  variety  of  models  are  available  for  the  description  of  the  laser-target  interaction  in  laser 
ablation.  In  this  study,  we  shall  model  laser  ablation  of  metals  in  a  vaeuum  using  ns-pulsed 
lasers  with  irradianees  in  the  10*- 104  W/em  range.  The  laser  irradianees  are  chosen  to 
produce  the  required  fluenee.  Several  processes  will  be  analyzed  in  this  model,  including 
surface  temperature  rise,  heat  conduction,  plasma  density,  ablation  rate  and  crater  depth. 


SURFACE  TEMPERATURE 


A  macroscopic  scale  description  of  the  laser-solid  interaction  is  used  in  this  study,  whereby 
the  thermal  heat  conduction  equation  is  used  to  analyze  heating  of  the  metal.  This  approach  is 
valid  for  laser  interaction  with  metals,  in  which  ease  laser  light  is  absorbed  by  interaction 
with  electrons,  and  is  then  transferred  to  lattice  phonons  by  collisions.  For  metals,  the  energy 
relaxation  time  is  of  the  order  of  10'n  s,  and  therefore  for  ns-pulsed  laser  beams  we  ean 
assume  optical  energy  is  turned  instantaneously  into  heat,  allowing  the  application  of  the  heat 
eonduetion  equation  \ 

The  temperature  rise  in  the  target  material  due  to  the  laser  beam  will  be  computed  from  the 
following  equation4, 

cp^af^=v^(r(2’°^r(2,/)^^w(2,0  (1) 

where  c  =  specific  heat  of  the  material 
p  =  material  density 

K  =  temperature  dependent  thermal  conductivity 
qiaser  —  heat  source  due  to  laser  radiation. 


The  direction  into  the  target  is  taken  as  the  positive-z  axis.  The  first  term  on  the  right-hand 
side  represents  the  heat  eonduetion  in  the  metal.  The  heat  source  term  in  Eq.  (1),  whieh  is  due 
to  the  absorption  of  laser  beam  radiation,  is  given  by 

<7w(z>0  =  /U(z)/(z,0  (2) 

where  p  =  absorption  coefficient 

I(z,t)  =  laser  beam  irradianee  as  a  function  of  time  and  position  in  the  target 


The  laser  beam  will  be  absorbed  by  the  metal  as  it  travels  through  the  target  material,  causing 
its  intensity  to  drop.  This  change  in  the  laser  beam  intensity  is  given  by  the  Beer-Lambert 
law4, 


(3) 


where  A  =  surfaee  absorptivity 

Io  =  irradianee  of  the  incident  beam 

Laser  irradianees  in  the  range  of  108-109  W/enT  will  be  considered  in  this  study.  The  effeet 
of  evaporation  of  the  metal  surfaee  is  incorporated  into  the  model  in  order  to  obtain  an 
aeeurate  description  of  the  thermal  field  in  the  target  material.  Conveetion  and  radiation 
effeets  will  be  included  in  the  future  to  better  describe  the  heating  effeet  of  the  laser. 


PLASMA  DENSITY 


The  plasma  density  at  the  surfaee  ean  be  determined  if  the  equilibrium  vapor  pressure  ean  be 
specified.  The  following  relation  gives  the  equilibrium  vapor  pressure  at  the  surfaee  as 

log,0  p  =  11.60  —  15880/  (4) 

/ 1  o  where  To  = 

temperature  at  the  surfaee 

The  equation  of  state  p=nokToean  then  be  used  to  ealeulate  the  density  at  the  surfaee. 


ABLATION  RATE 


The  ablation  rate  will  be  ealeulated  using  the  kinctie  model  described  in  Keidar  et  af.  A 
schematic  representation  of  the  plasma  and  kinetic  layers  near  the  surfaee  of  the  target  is 
shown  in  Figure  1.  In  the  kinctie  region,  using  Anisimov’s  assumption  that  the  velocity 
distribution  function  for  the  returned  partieles  is  Pfx(y),  where  f5  is  a  proportionality 

eocffieient,  the  relation  of  the  heavy  partiele  parameters  at  the  outer  boundary  of  the  kinctie 
layer  is  given  by  the  following  set  of  equations. 


n„ 


2(nd0) 


0-S 


=  n,p,  +  p 


2  (ndx) 


0-5 


X  {exp(— a2)  —  anoserfc(a)} 


"o 


4  d0  2  dx 


=  — --{(I  +  2a2)  —  /?[(0.5  +  a2)erfc(a )  —  a  exp(—a2)/nos]) 


(5) 


,  is  =  J  Vs  (a(g2  +  2-5)  “  0.5^[a(a2  +  2.5 )er/c(a)  -  (a:  +  2)  £>a73(-a-)/7r05]} 
(7rd0)ls  nd^ 

where  a  =  V,'1/(2kri/Tn)05 
=  m/(2kT0) 
dx  =  m/(2kTx) 
er/c(a)  =  1  —  er/(a) 
erf  (a)  =  error  function. 

For  relatively  small  plasma  densities,  we  ean  assume  that  the  plasma  accelerates  to  the  sound 
speed  at  the  edge  of  the  kinetie  layer,  and  in  this  ease  the  solution  of  the  kinetie  layer 


equations  (5)  gives  the  plasma  density  at  the  edge  of  this  layer  as  being  equal  to  0.34/7o  and 
the  temperature  as  0.67  To.  The  ablation  rate  is  then  given  by  the  relation1 1 

G  =  mnjl^  (6) 

where  m  =  mass  of  the  particles 
m  =  plasma  density 

Vi  =  velocity  at  the  edge  of  the  kinetie  layer 

Since  we  have  assumed  expansion  to  sound  speed,  the  velocity  at  the  edge  of  the  kinetic 
layer  is  then  given  by 

Vi  =  yjlkfjm 


n0  Hi 


Figure  1  Schematic  representation  of  kinetic  layer  near  the  surface 


LASER  PULSE 

The  laser  pulse  assumed  in  this  study  had  a  Gaussian  profile  with  FWHM  of  8  ns,  and 
centered  at  1 5  ns.  Furthermore,  a  KrF  laser  was  modeled,  with  wavelength,  X,  of  248  nm.  In 
this  study,  the  laser  fluencc  is  chosen  as  the  parameter,  and  the  laser  beam  intensity  is 
calculated  to  produce  the  required  fluencc.  Fluenees  in  the  range  of  3. 0-8.0  J/crn  were 
considered  in  this  paper.  An  example  of  the  laser  beam  profile  is  shown  in  Figure  2  for  the 
case  of  8.0  J/cm". 


Figure  2  Laser  beam  intensity  profile  used  in  the  model.  This  profile  corresponds  to  a 

fluence  of  8.0  J/cm2 


Table  1  Parameter  values  in  the  model  with  aluminum  as  the  target  material. 


Specific  Heat,  c 

940.0  J  kg'  K'1  (solid),  1289.0  J  kg'1  K'1 
(liquid) 

Mass  Density,  p 

2700.0  kg  m‘3  (solid),  2375.0  kg  m'7  (liquid) 

Melting  Point,  Tm 

933.5  K 

I  leat  of  Evaporation,  Hev 

10.8xl06  J  kg'1 

Atomic  mass,  m 

27  g  mol’1 

Absorption  Coefficient,  p 

5.7x1 07  m*' 

Absorptivity,  A 

0.21 

RESULTS  OF  THE  MODEL 

The  laser  model  was  applied  to  an  aluminum  target,  with  material  properties  as  given  in 
Table  1.  Aluminum  was  chosen  as  the  target  because  it  is  a  material  of  interest  and  to  allow 
for  comparison  of  the  results  with  other  models. 


TARGET  HEATING 

The  aluminum  target  is  initially  at  room  temperature  of  300  K.  As  a  result  of  the  incident 
laser  beam,  the  target  will  be  heated.  The  temperature  distribution  in  the  target  material  is 
plotted  at  several  representative  times  in  Figure  3.  As  expected,  the  temperature  is  maximum 
at  the  surface  for  all  times. 

Initially  the  temperature  near  the  surfaee  rises  quickly  as  the  laser  beam  intensity  reaehes  its 
highest  value  at  15  ns.  The  maximum  temperature  is  achieved  at  the  surfaee  several 
nanoseconds  after  the  peak  laser  intensity.  After  this  the  surfaee  temperature  as  well  as  the 
temperature  in  the  target  begins  to  drop  as  the  amount  of  energy  being  absorbed  is  reduced. 


The  temperature  gradients  in  the  target  also  decrease  and  a  smoother  profile  is  observable  at 
large  times.  It  can  also  be  observed  that  the  heat  conduction  has  penetrated  to  about  10  pm 
into  the  material  in  100  ns. 
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Figure  3  Calculated  temperature  profile  in  the  aluminum  target  for  a 
fluence  of  3.0  J/cm2  (a)  and  6.0  J/cm2  (b)  as  a  function  of  time 


SURFACE  TEMPERATURE 

The  surface  temperature  profile  as  a  function  of  time  is  plotted  in  Figure  4  for  several  laser 
fluenees.  Higher  laser  fluenee  results  in  a  higher  temperature  rise  at  the  surface,  as  expected. 
This  is  because  more  energy  is  absorbed  by  the  target  material.  Figure  4  (a)  shows  the 
surface  temperature  variation  with  time  for  fluenees  of  3.0  and  3.3  J/em:.  It  is  found  that  a 
maximum  temperature  of  2880  K  is  obtained  for  a  fluenee  of  3.3  J/enF  and  2560  K  for  the 
3.0  J/em  fluenee  ease.  After  reaching  a  maximum  about  6  ns  after  the  peak  laser  intensity, 
the  surface  temperature  drops  off  slowly  as  heat  is  transferred  into  the  aluminum  target. 


Similarly,  Figure  4  (b)  shows  that  maximum  surface  temperatures  of  3660  K,  5500  K  and 
6530  K  we  obtained  for  flucnccs  of  4.0,  6.0,  and  8.0  J/cm2  respectively. 

These  results  agree  with  results  in  Petcrlongo  ct  al.8,  where  a  triangular  laser  beam  profile 
was  used  instead.  A  maximum  surface  temperature  of  2588  K  is  obtained  about  5  ns  after  the 
peak  laser  intensity  for  a  fluencc  of  3.0  J/cnT.  For  3.3  J/em“,  the  maximum  surface 
temperature  obtained  was  2859  K.  Amoruso9  used  a  square  laser  beam  profile  to  study  the 
ablation  process  of  aluminum  and  obtained  surface  temperature  profiles  that  qualitatively 
looked  different  than  the  profiles  obtained  here.  This  was  due  to  the  different  laser  beam 
profiles  assumed.  Nevertheless,  maximum  surface  temperatures  of  about  2650  K  and  4400  K 
were  obtained  for  fluences  of  3.0  and  6.0  J/cnT.  Melc  ct  al.  meanwhile  used  a  KrF  laser 
with  Gaussian  profile  (FWHM  of  18  ns)  and  obtained  maximum  surface  temperatures  of 
about  3150  K  and  5200  K  for  fluences  of  3.0  and  6.0  J/em  respectively. 


Figure  4  Calculated  surface  temperature  profile  in  the  aluminum  target  for  a 
fluence  of  3.0  and  3.3  J/cm2  (a)  and  4.0,  6.0  and  8.0  J/cm2  (b) 


ABLATION  RATE 


The  heating  of  the  aluminum  target  will  cause  the  target  material  to  start  evaporating,  and  the 
ablation  rate  will  generally  increase  as  the  surface  temperature  increases.  The  ablation  rate  of 
the  aluminum  during  the  laser  irradianee  is  shown  in  Figure  A.  The  results  indicate  that  the 
ablation  rate  increases  sharply  as  the  target  surface  temperature  rises,  and  begins  dropping  as 
the  laser  beam  intensity  drops  and  the  surface  temperature  decreases.  In  fact,  the  vaporization 
of  material  is  one  of  the  faetors  eausing  the  surface  temperature  to  drop. 

The  results  for  ablation  rate,  shown  in  Figure  5  (a)  indicates  that  a  maximum  ablation  rate  of 
16.06  kg/rrTs  is  obtained  for  a  fluenee  of  3.0  J/errT,  while  an  ablation  rate  of  74.57  kg/m2s  is 
obtained  for  a  fluenee  of  3.3  J/em  .  Meanwhile,  Figure  5  (b)  shows  ablation  rates  of 
approximately  22800  and  59400  kg/m  s  are  obtained  for  fluenees  of  6.0  and  8.0  J/enr 
respectively.  Therefore  there  are  several  orders  of  magnitude  difference  in  the  ablation  rates 
for  fluenees  in  the  range  of  3.0  to  8.0  J/em  .  Nevertheless,  the  results  for  fluenees  above  6.0 
J/em  are  not  expected  to  be  accurate  as  both  AmorusoJ  and  Mele  et  al. 10  indicate  that  laser 
absorption  by  the  plasma  plume  becomes  significant  above  this  threshold. 

Peterlongo  et  al."  obtained  a  maximum  ablation  rate  of  18.9  kg/m‘s  for  a  fluenee  of 
3.0  J/em"  and  59.4  kg/m  s  for  a  fluenee  of  3.3  J/em\  This  difference  in  maximum  ablation 
rate  is  possibly  caused  by  the  different  methods  used  to  evaluate  the  rate  of  ablation.  The 
kinetic  approach  used  in  this  study  is  a  more  detailed  and  accurate  model  and  therefore 
should  give  a  better  approximation. 


Figure  5  Ablation  rate  as  a  function  of  time  for  fluenees  of  3.0  and  3.3  J/cm2  (a) 

and  4.0,  6.0,  8.0  J/cm2  (b) 


EVAPORATED  DEPTH 

The  laser  induced  ablation  of  the  aluminum  target  results  in  material  loss  from  the  surface.  In 
an  experimental  setting,  this  would  result  in  a  crater  around  the  laser  beam  incidence 
location.  In  this  1-D  study,  this  material  removal  causes  in  the  thickness  of  the  target  material 
to  be  reduced.  The  depth  of  evaporation  is  shown  in  Figure  6.  This  shows  most  of  the  change 


in  depth  occurring  between  12  ns  and  33  ns  after  the  beginning  of  the  laser  pulse  irradiation. 
After  that,  as  the  ablation  rate  drops,  the  evaporation  depth  levels  off. 

Fluenec  of  3.0  and  3.3  J/enr  produce  an  evaporation  depth  of  about  0.021  and  0.105  nm 
respectively,  as  shown  in  Figure  6  (a).  Figure  6  (b)  shows  that  fluenees  of  4.0,  6.0  and  8.0 
J/em2  produce  evaporation  depths  of  1.580,  45.00,  and  132.2  nm  respectively.  Again  there 
arc  several  orders  of  magnitude  difference  between  the  evaporation  depth  for  the  fluenees  in 
consideration 


Mele  ct  al.  0  obtained  evaporation  depths  of  about  6.5  and  48  nm  for  fluenees  of  4.0  and  6.0 
J/enV  respectively.  This  is  generally  in  agreement  with  results  of  this  present  model. 
Bogaerts  ct  al.  ‘  analyzed  laser  ablation  for  a  copper  target  and  found  that  for  laser  intensities 
in  the  range  of  10  -10  W/em  ,  the  evaporation  depth  was  in  the  range  of  0-100  nm. 


Figure  6  Evaporation  depth  as  a  function  of  time  for  fluenees  of  3.0  and  3.3  J/cm2  (a) 

and  4.0,  6.0,  8.0  J/cm2  (b) 


CONCLUSION 

This  report  has  described  a  model  for  nanosecond  pulsed  laser  ablation  of  an  aluminum  target 
in  vacuum  conditions  for  eases  when  the  laser  beam  absorption  by  the  plasma  above  the 
target  is  not  significant.  This  model  describes  the  laser-solid  interaction  that  results  in  target 
heating  and  vaporization  leading  to  ablation.  The  results  include  the  temperature  distribution 
in  the  aluminum  target,  the  surface  temperature  profde,  the  ablation  rate  and  the  amount  of 
evaporation. 

Laser  fluenees  in  the  range  of  3. 0-8.0  J/em  were  considered.  Studies  indicate  that  for 
fluenees  above  this  range,  laser  absorption  by  the  plasma  created  by  the  ablation  becomes 
significant  and  has  to  be  taken  into  account.  Results  show  that  temperatures  of  2560  K  and 
5500  K  are  obtained  for  fluenees  of  3.0  and  6.0  J/em  .  Furthermore,  0.021  nm  of  aluminum 
was  evaporated  for  a  fluenec  of  3.0  J/enT  while  45  nm  was  evaporated  for  a  fluenee  of  6.0 
J/em  .  The  results  show  acceptable  agreement  with  results  of  other  analytical  ablation 


models,  although  the  use  of  different  laser  beam  profiles  makes  any  accurate  comparison 
difficult. 


It  should  be  reiterated  that  only  after  considering  the  plasma  plume  expansion  above  the 
target  can  accurate  results  be  obtained  for  the  higher  fluenees  whieh  arc  of  interest. 
Nevertheless,  this  goal  of  this  report  was  to  demonstrate  that  the  laser-solid  interaction  model 
presented  here  is  acceptable  for  lower  fluenees.  The  results  obtained  using  this  ablation 
model,  particularly  the  temperature,  velocity  and  density  at  the  edge  of  the  kinctie  layer  will 
be  used  as  boundary  conditions  to  model  the  plasma  plume  expansion  above  the  target  in 
order  to  obtain  accurate  results  for  laser  ablation  at  higher  fluenees. 


REFERENCES 

Radzicmski,  L.J.  and  Cremcrs,  D.A.,  Laser-Induced  Plasmas  and  Applications.  Marcel 
Dekker,  Inc.  (1989) 

:Keidar,  M.,  Boyd,  l.D.  and  Beilis,  1. 1.,  “On  the  model  of  Teflon  ablation  in  an  ablation 
controlled  discharge.”  J.  Phvs.  D:  Appl.  Phvs.  34.  (2001)  1675 

'von  Allmen,  M.,  Laser  Beam  Interactions  with  Materials.  Springer  Heidelberg  (1987) 

1  Semak,  V.V.  and  Miller,  T.F.,  “Modeling  of  Laser  Charring  and  Material  Removal  in 
Fiberglass  Materials,”  J.  Directed  Energy  2(1),  (2006)  5 
Dushman,  S.,  Vaeuum  Techniques,  John  Wiley  &  Sons,  Inc.  (1949) 

6  Keidar,  M.,  Fan,  J.,  Boyd,  l.D.  and  Beilis,  I.I.,  “Vaporization  of  heated  materials  into 
discharge  plasmas,”  J  Appl.  Phvs.  89(6),  (2001)  3095 
Anisimov.  S.l.  Sov.  Phvs.  JETP  27.  (1968)  182 
'  Pctcrlongo,  A.  ,  Miotcllo,  A.,  and  Kelly,  R.,  “Laser-pulse  sputtering  of  aluminum: 
Vaporization,  boiling,  superheating,  and  gas-dynamic  effeets,”  Phvs.  Rev.  E  50(6),  (1994) 
4716 

Amoruso,  S.,  “Modeling  of  UV  pulsed-lascr  ablation  of  metallic  targets,”  Appl.  Phvs.  A  69, 
(1999)323 

10  Mele,  A.,  Guidoni,  A.G.,  Kelly,  R.,  Flamini,  C.  and  Orlando,  S.,  “Laser  ablation  of  metals: 
Analysis  of  surface-heating  and  plume-expansion  experiments,”  Appl.  Surf.  Sci.  109,  (1997) 
584 

nYilbas,  B.S.,  “Numerical  approach  to  pulsed  laser  heating  of  Semi-Infinite  Aluminum 
substance,”  Heat  and  Mass  Tran_s.  31,  (1996)  279 

'  "Bogacrts,  A.,  Chen,  Z.,  Gijbcls,  R.  and  Vertes,  A.,  “Laser  ablation  for  analytical  sampling 
what  ean  we  learn  from  modeling?,”  Spectrochimica  Aeta  Part  B  58,  (2003)  1867 
nKeidar,  M.,  Boyd,  l.D.,  Luke,  J.  and  Phipps,  C.,  “Plasma  generation  and  plume  expansion 
for  a  transmission-mode  microlaser  ablation  plasma  thruster,”  J.  Appl.  Phvs.  96(1),  (2004)  49 
i4Boardman,  A.D.,  Crcsswcll,  B.  and  Anderson,  J.,_“An  analytical  model  for  the  laser 
ablation  of  materials,”  Appl.  Surf.  Sei.  96-98,  (1996)  55 

|SLebocuf,  J.N.,  Chen,  K.R.,  Donato,  J.M.,  Geohcgan,  D.B.,  Liu,  C.L.,  Purctzky,  A. A.  and 
Wood,  R.F.,  “Modeling  of  dynamical  processes  in  laser  ablation,”  Appl.  Surf.  Sei.  96-98, 
(1996)  14 

''Tannchill,  J.C.,  Anderson,  D.A.  and  Pletehcr,  R.H.,  Computational  Fluid  Mechanics  and 
Heat  Transfer,  Taylor  &  Francis  (1997) 


